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Abstract 

A random field composed by Poisson distributed Brownian vortex filaments is 
constructed. The filament have a random thickness, length and intensity, governed 
by a measure 7. Under appropriate assumptions on 7 we compute the scaling law 
of the structure function and get the multifractal scaling as a particular case. 

1 Introduction 

Isotropic homogeneous turbulence is phenomenologically described by several theories, 
which usually give us the scaling properties of moments of velocity increments. If u{x) 
denotes the velocity field of the fluid and Sp (e), the so called structure function, denotes 
the p-moment of the velocity increment over a distance e (often only its longitudinal 
projection is considered), then one expect a behavior of the form 

Spie) = {\uix + e)-uix)\P) ^e'^^. (1) 

Here, in our notations, e is not the dissipation energy, but just the spatial scale parameter 
(see remarking. Let us recall two major theories: the Kolmogorov-Obukov scaling law 
(K41) (see [H]) says that 

probably the best result compared with experiments; however the heuristic basis of the 
theory also implies = f which is not in accordance with experiments. Intermittency 
corrections seem to be important for larger p's. A general theory which takes them into 
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account is the multifractal scaling theory of Parisi and Frisch [llj, that gives us (p in the 
form of a Fenchel-Legendre transform: 

Cp = mf[hp + 3- D{h)]. 

This theory is a sort of container, which includes for instance the striking particular case 
of She and Leveque jlHj . We do not pretend to go further in the explanation of this topic 
and address the reader to the monograph [Hj. 

The foundations of these theories, in particular of the multifractal one, are usually 
mathematically poor, based mainly on very good intuition and a suitable "mental image" 
(see the beginning of Chapter 7 of P). Essentially, the scaling properties of Sp{e) are 
given a priori, after an intuitive description of the mental image. The velocity field of the 
fluid is not mathematically described or constructed, but some crucial aspects of it are 
described only in plain words, and then Sp{e) is given (or heuristically "deduced"). 

We do not pretend to remedy here to this extremely difficult problem, which ultimately 
should start from a Navier-Stokes type model and the analysis of its invariant measures. 

The contribution of this paper is only to construct rigorously a random velocity field 
which has two interesting properties: i) its realizations have a geometry inspired by the 
pictures obtained by numerical simulations of turbulent fiuids; ii) the asymptotic as e —>■ 
of Sp{e) can be explicitly computed and the multifractal model is recovered with a suitable 
choice of the measures defining the random field. Its relation with Navier-Stokes models 
and their invariant measures is obscure as well (a part from some vague conjectures, see 
[Zj), so it is just one small step beyond pure phenomenology of turbulence. 

Concerning (i), the geometry of the field is that of a collection of vortex filaments, 
as observed for instance by 19] and many others. The main proposal to model vortex 
filaments by paths of stochastic processes came from A. Chorin, who made several con- 
siderations about their statistical mechanics, see The processes considered in j3] are 
self-avoiding walks, hence discrete. Continuous processes like Brownian motion, geomet- 
rically more natural, have been considered by [THI, [TE], |E], [H], [HI and others. We do 
not report here numerical results, but we have observed in simple simulations that the 
vortex filaments of the present paper, with the tubular smoothening due to the parameter 
i (see below), have a shape that reminds very strongly the simulations of [3^. 

Concerning (ii), we use stochastic analysis, properties of stochastic integrals and ideas 
related to the theory of the Brownian sausage and occupation measure. We do not know 
whether it is possible to reach so strict estimates on Sp{e) as those proved here in the 
case when stochastic processes are replaced by smooth curves. The power of stochastic 
calculus seem to be important. 

Ensembles of vortex structures with more stiff or artificial geometry have been consid- 
ered recently by PP and ^2] • They do not stress the relation with multifractal models and 
part of their results are numerical, but nevertheless they indicate that scaling laws can be 
obtained by models based on many vortex structures. Probably a closer investigation of 
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simpler geometrical models like that ones, in spite of the less appealing geometry of the 
objects, will be important to understand more about this approach to turbulence theory. 
Let us also say that these authors introduced that models also for numerical purposes, so 
the simplicity of the structures has other important motivations. 

Finally, let us remark about a difference with respect to the idea presented in jSj and 
also in There one put the attention on a single vortex and try to relate statistical 
properties of the path of the process, like Flory exponents of 3D self-avoiding walk, with 
scaling law of the velocity field. Such an attempt is more intrinsic, in that it hopes to 
associate turbulent scalings with relevant exponents known for processes. On the contrary, 
here (and in and f^) we consider a fluid made of a multitude of vortex structures and 
extract statistics from the collective behavior. In fact, it this first work on the subject, 
we consider independent vortex structures only, having in mind the Gibbs couplings of 
[S] as a second future step. Due to the independence, at the end again it is just the single 
filament that determines, through the statistics of its parameters, the properties of Sp{e). 
However, the interpretation of the results and the conditions on the parameters are more 
in the spirit of the classical ideas of K41 and its variants, where one thinks to the 3D 
space more or less filled in by eddies or other structures. It is less natural to interpret 
multifractality, for instance, on a single filament (although certain numerical simulation 
on the evolution of a single filament suggest that multifractality could arise on a single 
filament by a non-uniform procedure of stretching and folding [Sj). 

1.1 Preliminary remarks on a single filament 

The rigorous definitions will be given in section 2. Here we introduce less formally a few 
objects related to a single vortex filament. 

We consider a 3d-Brownian motion {Xt}te[o,T] starting from a point Xq. This is the 
backbone of the vortex filament whose vorticity field is given by 



The letter t, that sometimes we shall also call time, is not physical time but just the 
parameter of the curve. All our random fields are time independent, in the spirit of 
equilibrium statistical mechanics. We assume that pe{x) = p{x/i) for a radially symmetric 
measurable bounded function p with compact support in the ball -5(0, 1) (the unit ball in 
3d Euclidean space). To have an idea, consider the case p = 1b{o,i)- Then ^single (a^) = 
outside an £-neighboor Ui of the support of the curve {Xt}tg[o,r]- Inside Ue, ^singic(a;) is a 
time-average of the "directions" dXt, with the pre-factor U/i'^. More precisely, if a; e Ui, 
one has to consider the time-set where X^ e B{x,i) and average dXf on such time set. 
The resulting field ^singie(a^) looks much less irregular than {X^j^g^.T], with increasing 
irregularity for smaller values of i. 
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When p is only measurable, it is not a priori clear that the Stratonovich integral 
^singie(a;) is well defined, since the quadratic variation corrector involves distributional 
derivatives of p (the Ito integral is more easily defined, but it is not the natural object to 
be considered, see the remarks in (HI). Since we shall never use explicitly ^single (a;), this 
question is secondary and we may consider ^single (a^) just as a formal expression that we 
introduce to motivate the subsequent definition of 'Usingio(a;). However, at least in some 
particular case (p = 1b(o,i)) ^r under a little additional assumptions, the Stratonovich 
integral (^single (a;) is well defined since the corrector has a meaning (in the case of p = 1 b(o,i) 
the corrector involves the local time of 3d Brownian motion on sferical surfaces). 

The factor U /^"^ in the definition of ^singie(a^) is obscure at this level. Formally, it could 
be more natural just to introduce a parameter F, in place of t//£^, to describe the intensity 
of the vortex. However, we do not have a clear interpretation of F, a posteriori, from out 
theorems, while on the contrary it will arise that the parameter U has the meaning of 
a typical velocity intensity in the most active region of the filament. Thus the choice of 
the expression U / has been devised a posteriori. The final interpretation of the three 
parameters is that T is the length of the filament, £ the thickness, U the typical velocity 
around the core. 

The velocity field u generated by ^ is given by the Biot-Savart relation 

U /"^ 

Msmgic(a;) = ^ j ^ii^ - ^t) A odXt (3) 
where the vector kernel Ki{x) is defined as 

= Vl'^W = ^ / P,{y)^,iy, Vd^) = ~[ (4) 

47r JB(o,e) \x - 47r ^^(o,^) \x - y\ 

The scalar field Ve satisfy the Poisson equation AVe = pi in all M^. Since p is radially 
symmetric and with compact support we can have two different situations according to 
the fact that the integral of p: Q = J^^q ^-j p{x)dx is zero or not. If Q = then the field 
is identically zero outside the ball 5(0, £). Otherwise the fields V^, i^'^ outside the ball 
-8(0, £) have the form 



X 



V,{x) = Q^, Ke{x) = -2Qef-, \x\ > L (5) 



\x\ \x 



Accordingly we will call the case Q = Q short range and Q 7^ long range. The proof 
of the previous formula for K is given, for completeness, in the remark at the end of the 
section. 

A basic result is that the Stratonovich integral in the expression of Usingie can be 
replaced by an Ito integral: 
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Lemma 1 ltd and Stratonovich integrals in the definition of Usmgiei^) coincide: 



where the integral is understood in ltd sense. Then Usmgie is a (local) -martingale with 
respect to the standard filtration of X . 

About the proof, by an approximation procedure that we omit we may assume p 
Holder continuous, so the derivatives of Ki exist and are Holder continuous by classical 
regularity theorems for elliptic equations. Under this regularity one may compute the 
corrector and prove that it is equal to zero, so, a posteriori, equation ^ holds true in the 
limit also for less regular p. About the proof that the corrector is zero, it can be done 
component-wise, but it is more illuminating to write the following heuristic computation: 
the corrector is formally given by 



Since is a gradient, we have curl = 0, so the corrector is equal to zero. 

Remark 1 One may verify that div Ugmgie = 0, so Ugmgie "i^dy be the velocity fluid of an 
incompressible fluid. On the contrary, div ^single is different from zero and curl Usingie is 
not C,singie but its projcction on divergence free fields. Therefore, one should think of single 
as an auxiliary field we start from in the construction of the model. 

Remark 2 Let us prove limited to K to avoid repetitions. We want to solve /W = p 
with p spherically symmetric. The gradient K = W of V satisfies V ■ K = p and, by 
spherical symmetry, it must be such that 






Now, from the property dX^dX^ = 6ijdt one can verify that 



{VKe{x - Xt)dXt) A dXt = (curl K^) {x - Xt)dt. 



K{x) 



X 



X 



for some scalar function f{r). By Gauss theorem we have 




where da{x) is the outward surface element of the sphere. So 
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namely f{r) = Therefore 

K{x)^Q{r) 



47r|xP 

with Q{r) — Q{1) if r > 1 (since p has support in B{0, 1)). 

2 Poisson field of vortices 

Intuitively, we want to describe a collection of infinitely many independent Brownian vor- 
tex filaments, uniformly distributed in space, with intensity-thickness-length parameters 
([/, £, T) distributed according to a measure 7. The total vorticity of the fluid is the sum 
of the vorticity of the single filaments, so, by linearity of the relation vorticity-velocity, 
the total velocity field will be the sum of the velocity fields of the single filaments. 
The rigorous description requires some care, so we split it into a number of steps. 

2.1 Underlying Poisson random measure 

Let 5 be the metric space 

H = {(?7, £, T, X) e X C([0, 1]; M^) : < £ < Vt < 1} 

with its Borcl cr- field E (S). Let (^2, P) be a probability space, with expectation denoted 
by and let fi^, a; G f2, be a Poisson random measure on B (S), with intensity v (a a- 
finite measure on (S)) given by 

dv{U,l,T,X) ^ d^{U,l,T)d'W{X). 

for 7 a (T-finite measure on the Borel sets of {(t/, T) e : < ^ < < 1} and 
dyV{X) the cr-finite measure defined by 



/ i:{X)dW{X) = 

Jc([0,l];K3) 



ij{x)m,,{x) 

C([0,l];K3) 



dxQ 



for any integrable test function : C([0, 1]; R^) R where dW^oi^) is the Wiener mea- 
sure on C([0, 1],R^) starting at Xq and dx^ is the Lebesgue measure on M?. Heuristically 
the measure W describe a Brownian path starting from an uniformly distributed point in 
all space. The assumptions on 7 will be specified at due time. 

The random measure //^^ is uniquely determined by its characteristic function 

Eexp [i j^viiMdi)^ = exp (-^(^^^^^ - 1)^(^0 
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for any bounded measurable function 9? on S with support in a set of finite i/-measure. 
In particular, for example, the first two moments of /i read 



E 



and 



We have to deal also with moments of order p; some useful formulae will be now given. 
2.2 Moments of the Poisson Random Field 

Let : 5 ^ M be a measurable function. We shall say it is /^-integrable if it is |J,^^- 
integrable for P-a.e. a; e In such a case, by approximation by bounded measurable 
compact support functions, one can show that the mapping a; h- > /^^ ((^) is measurable. 

If yj' : S M is a measurable function with z/ ((/? 7^ (/?') = 0, since P {^{ip ^ (/?') = 0) = 
1, we have that is /i-integrable and P{^{ip) = /i ((/?')) = 1. Therefore the concept of 
/x-integrability and the random variable 11 (ip) depend only on the equivalence class of ip. 

Let ip he a measurable function on H, possibly defined only u-a..s. We shall say that it 
is ^-integrable if some of its measurable extensions to the whole S is //-integrable. Neither 
the condition of /x-integrability nor the equivalence class of {ip) depend on the extension, 
by the previous observations. We need all these remarks in the sequel when we deal with 
ip given by stochastic integrals. 

Lemma 2 Let ip he a measurable function on S (possibly defined only v-a.s.), such that 
^{(fi^) < 00 for some even integer number p. Then (p is /i-integrable, ii{(f) e (D,), and 



Proof. Assume for a moment that ip is bounded measurable and with support in a 
set of finite //-measure, so that the quahtative parts of the statement are obviously true. 
Using the moment generating function 



If in addition we have ^{ip^) — for every odd k < p, then 



E[^{ipf]>v{^^). 
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we obtain 



OO ,p oo ^ 

p=0 ^' n=0 



oo 



n=l 

oo 



A 



k=l 



n\ ^ — ' ki\ ■ ■ ■ kj. 

n=l fci,...,fc,,>l 



CO p 
p=l ^ n=l A:i,...,A;„>l ^ " 

fclH |-fcn=P 

Hence we have an equation for the moments: 

p 



n=l fci,...,A;„>l 
fciH hfcn=P 



nl/ci! ■ ■ ■ kj. 



Since i/dy^l'^) < for k < p we have 

E[iM^)n<E[Mi^in 



< 



E E 

n=l A;i,...,fe„>l 
fciH \-kn=p 



p\ 



n\ki\ ■ ■ ■ kr 



p\ 



^''(M')E E „,,,,..,,, 

n=l fci,...,fc„>0 
fciH |-A:ti=P 

P „ 

— ^ n! 

n=l 

This proves the first inequahty of the lemma. A posteriori, we may use it to prove the 
quahtative part of the first statement, by a simple approximation procedure for general 
measurable ip. 

For the lower bound, from the assumption that = for odd < p, in the sum 

((Zj) we have contributions only when all ki, i = 1, . . . ,n are even. Then, neglecting many 
terms, we have 

E[/i(^)^] > ui^n- 

The proof is complete. 
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2.3 Velocity field 

Let p be a radially symmetric measurable bounded function p on with compact support 
in the ball -8(0, 1), and let Ki be defined as in section ITTTl Then we have that 

Ki is Lipschitz continuous. (8) 

Indeed by an explicit computation it is possible to show that 

K,{x) = -2Q{\x\) 



\x 



3 



with 



Q{r) := / p{x)dx 

' B(Q,r) 



(since p has support in 5(0, 1) we have Q{r-) = Q if r > 1). Then 



VK,{x) = -2Q'{\x\)^-2Q{\x\) 



\x 



1 X® X 



with Q'{r) = dQ{r)/dr and we can bound 

\VKi{x)\<C\\p\\^ 

since 

Q(r) < C||p|Ur^ Q'ir)<C\\p\\^r' 



For any x,Xq G M^, i,T > 0, the random variable X t—* K£{x — Xt) A dXt is 
defined W^^-a-s. on ^([0,1],^^). We also have, given a; G M^ £,T > 0, that X ^ 
J^' Ke{x-Xt)AdXt is a well defined measurable function, defined W-a.s. on C([0, 1], M'^). 
More globally, writing ^ = {U, i, T, X) for shortness, for any x G we may consider the 
measurable M^-valued function 



U 



T 



e ^ 4ngle (^) := Y2 ^^^(^ - ^ 

defined z/-a.s. on S. In plain words, this is the velocity field at point x of a filament 
specified by ^. 

Again in plain words, given uj E fl, the point measure p^^ specifies the parameters and 
locations of infinitely many filaments: formally 



p = J2^^- (9) 
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for a sequence of i.i.d. random points {^°'} distributed in S according to u (this fact is 
not rigorous since u is only a-finite, but it has a rigorous version by localization explained 
below). Since the total velocity at a given point x G should be the sum of the 
contributions from each single filament, i.e. in heuristic terms 



U{X) = 



E«Singlc(^) (10) 

we see that, in the rigorous language of /i, we should write 

u{x) = / 4^gi^ (x) fi{d^) = (w^i^gig (x)) . (11) 



If we show that u (|Msingic < oo for some even p > 2, then from lemma ^single (^) 

is yU-integrable and the random variable 



uj ^ u{x, uj) := /i^ {u,ingie {x)) 



is well defined. 

Since the proof of the next lemma we start to use the occupation measure of the 3d 
Brownian motion. We set 





for every Borel set B of M^. 

Lemma 3 Given x G and p > 0, there exist Cp > such that for every i"^ < T < Iwe 
have 

W [W{xY^^] < CpfPiT 

where 

Wix)= / \Keix - Xt)\'' dt. 
Jo 

Proof. Let us bound K hj a multi-scale argument. This is necessary only in the long- 
range case (see the introduction). If \y\ < i we can bound \Ki{y)\ < C£. Indeed if \w\ < 1, 
1/^(^)1 < C for some constant C and then if \y\ < i we have \Ki{y) \ = i\Ki{y/i)\ < Ci. 

Next, given A > i and an integer A^, consider a sequence {^i}j=o '^^ scales, with 
i = io < £i < ... < £n = A. Then, for z = 1, . . . , A^, if ii^i < \y\ < ii, by the explicit 
formula for Ki{y) we have \Ki{y)\ < Ci{i/\y\y. Therefore 

my)\<ce(^-^^ Wm<.. 

If \y\ > A we simply bound K^{y) < Ci{i/A)'^. 
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Summing up, 

\Keiy)f = \Keiy)f [ l\y\<e + Vi<|y|<^i + 1m>a 



N 



i=l 



i=i V » 1 
which imphes the following bound for W{x): 

N 



hi-i<\y\<ei + Ci [-] l|y|>A 



i=l 



where has been defined above. By the additivity of S i-^ L^, the sum appearing in 
this equation can be rewritten as 

N / n \ i N-1 



E 

1=1 



ii-i 



L'B{xA)\Bix,ei-i) ~ X] 
i=l 



+ 



4 



^B{x,ii) 



£^_^ I ~B{x,A) ^B{x, 



N-l 

1=1 



'-i-i 



^B{xA) + 



^B{x,A)- 



Assume that ii/ii^i < 2 uniformly in i = 1, . . . , A^. Then 



i-l 



N 

E 

1=1 

N-l 

s^E 

i=l 

Notice now that 



J^B{xA)\B{xA-i) 



U-1 



^1 ^Wo + ^l^l ^ 



N-l 

E 

1=1 



'-i-i 



A 



J ^B{x,A)- 



< 1 



SO that, by Cauchy-Schwartz and Jensen inequalities we have 

lP/2 



J—^J ^B{x,£i)\B{x,ii-i) 



E 



i=l 

N-l 

i=l 



U-1 



p/2 
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An upper bound for W{xy^'^ is then obtained as 



p/2 



N-1 



i=l 



< 



A 



2p 



i=l 



'H-l 



where we have used again Cauchy-Schwartz inequahty. 

We use now lemma IT^ with a = 1. For a given A G (0, 1), we take bot e and ^ equal 
to A in (j21I) and ((23), and get 



Then we obtain 



< a 



2p 



N-l 



i=l 



and taking the limit as the partition gets finer: 



< Cpf nT + CpF ( - J 



2p 



u 



u 
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The integral can then be computed as 



/ - iu^VTYiuy Vffu^du 

Je JVt 



= FfTir^ -Vf '] + (p - i)-irf/'£f+2[r(i-p)/2 _ A(i-f)] 

< enr + (p - i)-^ep-^iT 
Vf 

Using the fact that £ < VT and letting A — > oo we finally obtain the claim. The proof is 
complete. □ 



Remcirk 3 The multiscale argument above can be rewritten in contunuum variables from 
the very beggining by means of the following identity: if f : [0, oo) W is of class and 
has a suitable decay at infinity, then 



/ f{\x-X,\)dt = - f (r) Ll^^^^^dr. 
Jo Jo 



This identity can be applied to W{x). The proof along these lines is not essentially shorter 
and perhaps it is more obscure, thus we have choosen the discrete multiscale argument 
which has a neat geometrical interpretation. 

Corollary 1 Assume 

-fiUPET) < oo 

for some even integer p > 2. Then, for any x eM.^, we have 

{\Kingle{x)n < ^ 

and the random variable 
has finite p-moment: 

E[\u{x)\''\<oo. 

Proof. We have 

^ (hsmgle(a;)r) = J ^xoiHingleix)\^)d7iU,i,T)dXo. 
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By Burkholder-Davis-Gundy inequality, there is Cp > such that 

TJP 



Hence 



jjp 



HP 

< C'p—fPiT = C'^UHT. 

Therefore v (|ttsingie i^)^) < oo by the assumption 'y{U^iT) < oo. The other claims are a 
consequence of lemma |21 □ 

Lemma 4 Under the previous assumptions, the law of u{x, ■) is independent of x and is 
invariant also under rotations: 

Ru{x, ■) = u{Rx, ■) 

for every rotation matrix R. 

Proof. With the usual notation ^ = ([/, i, T, X) we have 



single 



i^) = '^^ i^) = u'i'nS'^ (0) = "Idle (0) 



where Tx{U,i,T,X) = {U,i,T,X — x). The map is a measurable transformation of S 
into itself. One can see that u is Tx-invariant; we omit the details, but we just notice that 
u is not a finite measure, so the invariance means 



for every ip ^ (S, u). From this invariance it follows that the law of the random measure 
H is the same as the law of the random measure r^/i. Therefore 



■"single (0) 



^xfJ^) 4ngle(0) 



single 



(0) 



This proves the first claim. 

If is a rotation, from the explicit form of Ki it is easy to see that 



RK^iy) = K,{Ry) 
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hence 



^"Single i^) 



^ J RK^{x-Xt)AdRXt 

T 



u 



I KeiRx - RXt) A dRXt 
Jo 



u 



single 



{Rx) 



where we have set R{U,i,T,X) = {U,i,T, RX). Again Ru = u, Rfi = fi, so the end of 
the proof is the same as above. □ 
We say that a random field u{x,-) is homogeneous if it is independent of x and it is 
isotropic if it is invariant under rotations. 



-f{UP£T) < oo 

is an isotropic homogeneous random field, with 



Corollary 2 Assume 

for every p > 1. Then {u{x, G 
finite moments of all orders. 

This corollary is sufficient to introduce the structure function and state the main results 
of this paper. However, it is natural to ask whether the random field {u{x, ■);x G M.^} 
has a continuous modification. Having in mind Kolmogorov regularity theorem, we need 
good estimates of E [\u{x) — u{y)f]. They are as difficult as the careful estimates we shall 
perform in the next section to understand the scaling of the structure function. Therefore 
we anticipate the result without proof. It is a direct consequence of Theorem ^ 

Proposition 1 Assume 

7([/P£T) < oo 

for every p > 1. Then, for every even integer p there is a constant Cp > such that 



E [\uix) - uiy)f] < 



' A\x — y\ 



iT 



Consequently, if the measure 7 has the property that for some even integer p and real 



number a > 3 there is a constant Cp > such that 



7 



< C'e'' for any e G (0, 1) 



(12) 



then the random field u{x) has a continuous modification. 
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Remark 4 A sufficient condition for M^) is: there are a > 3 and /5 > such that for 
every sufficiently large even integer p there is a constant Cp > such that 

7 [Unr ■ le<ei-p] < Cpe"" for any e G (0, 1) . 



Indeed, 



1 



jjp 



iAe 



iT 



1 



JJP 



IT -1 



+ 7 



JJP 



so we have M^) with a suitable choice of p. This happens in particular in the multifractal 
example of section VJ/<^ remark\^ 

The model presented here has a further symmetry which is not physically correct. This 
symmetry, described in the next lemma, implies that the odd moments of the longitudinal 
structure function vanish, in contradiction both with experiments and certain rigorous 
results derived from the Navier-Stokes equation (see jH]). The same drawback is present 
in other statistical models of vortex structures fT]. 

Beyond the rigorous formulation, the following property says that the random field 
ig has the same law as —u'^i^ ^^- We cannot use the concept of law since ^ does not 



u 



sing. 



live on a probability space. 

Lemma 5 Given xi, ...,Xn G IR'^, the measurable vector 



U 



single 



has the property 



for every f = (f . : 1^'^" — * K with a polynomial bound in its variables.. 



Proof. With the notation X+ 



-u, 



single 



{Xk) 



Xx-ti we have 
U 



X,) A odX, 



u_ 



Ke{xk - X,) A odXs 



where S{U, £,T, X) = {U,i,T,X). Since Su = u, we have the result (using the integra- 
bility of corollary P). □ 
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Lemma 6 If p is an odd positive integer, then 

for every x,y eM.^. 

Proof. It is sufficient to apply the lemma to the function 

LP {ui, U2) := {ui -U2,y- xf 
and the points Xi = y, X2 = x, with the observation that 

ip {-ui, -U2) = -ip {Ui, U2) . 

□ 

2.4 Localization 

At the technical level, we do not need to localize the cr-finite measures of the present work. 
However, we give a few remarks on localization to help the intuitive interpretation of the 
model. Essentially we are going to introduce rigorous analogs of the heuristic expressions 
dni) and (fTn|) written at the beginning of the previous section. The problem there was that 
the law of should be z/, which is only a a-finite measure. For this reason one has to 
localize v and ^. 

Given A E B (H) with < z/ (A) < 00, consider the measure /xa defined as the 
restriction of fi to A: 

/iA (B) =fi{AnB) 

for any B G -6(2). It can be written (the equality is in law, or a.s. over a possibly 
enlarged probability space) as the sum of independent random atoms each distributed 
according to the probability measure B E B (A) ua {B) := z/(i?|y4): 

Na 

where Na is a Poisson random variable with intensity i^iA) and the family of random 
variables {^°}aeN is independent and identically distributed according to ua- Moreover if 
{AjjjgN is a family of mutually disjoint subsets of E then the r.v.s {/iAijieN are indepen- 
dent. 

Sets A as above with a physical significance are the following ones. Given < r/ < 1 
and R> Q, let 

A,,R = {iU,TJ,X)eE:i>T], \xo\ < R}. 
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In a fluid model we meet these sets if we consider only vortexes up to some scale rj (it 
could be the Kolmogorov dissipation scale) and roughly confined in a ball of radius R. 
If we assume that the measure 7 satisfies < ^{i > rj) < 00 for each 77 > 0, then 
< v (Ajj^r) < 00, and for the measure iiri,R '■— /J'A^ ^ we have the representation 

a=l 

where Nr,,R is V (i/(>l^,ij)) and {C"}aeN are i.i.d. with law Ur,,R ^a^.a- 

For any x e we may consider ^ 1— > wfj^gig (x) as a random variable in M^, defined 
Ujj^R-a.s. on Arj^R. Moreover we may consider the r.v. Ur,,R{x) on {fl, A, P) defined as 

Ur^A^)-= I 4ngle (^) /^(^O = / 4ngle (^) /^r,,il(c?0 • (13) 

It is the velocity field at point x, generated by the vortex filaments in Aj^^r e B (5). In 
this case we have the representation 

a=l a=l > 

where the quadruples = ([/°, T", X") are distributed according to u^^r and are 
independent. If a; i— > ^ (cu) is any one of such quadruples, the random variable 

^ ^ '"single (^) 

is well defined, since the law of ^ is v^^^r and the random variable ^ i— > Mangle (2^) is well 
defined z/^^^-a-S. on A^^ r. Therefore u^^A^) ^ ^^^^ defined random variable on (Q, A, P). 
We have noticed this in contrast to the fact that the definition of u{x) required difficult 
estimates, because of the contribution of infinitely many filaments. 

Given the Poisson random field, by localization wc have constructed the velocity fields 
Uri^R^x) that have a reasonable intuitive interpretation. Connections between Ur,^R{x) and 
u{x) can be estabhshed rigorously at various levels. We hmit ourselves to the following 
example of statement. 

Lemma 7 Assume 

-fiU^iT) < 00 

for any p > 0. Then, at any a; e R^, 

lim E[\ur,R(x) -u(xW]^0. 
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Proof. Since 



Un,R{x) - u{x) 



we have 



- 1) 



^single 



Notice that we do not have v (Aj^ ^) — 0, in general, so the argument to prove the lemma 
must take into account the properties of the r.v. Msingie(^)- We have (by Burkholder- 
Davis- Gundy inequality) 



"^single l-^/ 



ri 



^single (-^ 

\xo\>R Jl>ri 

J\xo\>R Jl>ri ^ 



f l^d^ ^' ^) / [Wixf/'] dxo 

J J\xo\>R 

where W{x) has been defined in lemma 01 Let us show that 

/ >v^o [wixy/^] dxo < Cpfnr ■ e (r) 

J\xo\>R 



(14) 



where 6 (R) as R ^ oo. The proof will be complete after this result. Recall from 
the proof of lemma El that we have 



2p 



N-1 



^From lemma 1131 we have 



i=l 



'i-l 



i^'B{x/i)T^'^- 



\xo\>R 



dxo 



(t'^ V 

< Cp(A A VTY\{\ V Vff exp ^_ -^-(^+^) 
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Hence 



\xo\ >R 
2; 



< CptnTexp 



[Wixf^^] dxo 
R-(\x\ + 



N-l 



1=1 



T 

2 / \ 2' 



A 



2p 



rpp/2 



R-{\x\ 



T 



\xo\ >R 



Repeating the arguments of lemma El we arrive at 

>v.„ [w{xy/^] dxo 

Since T and £ are smaller than one, and p > 2, we also have 

[Wixf/^] dxo 



'\xo\ >R 

< CpfPiT (exp {-R - (|x| + A)) + A-^p) 
With the choice A = R/2 we prove (fT^ . The proof is complete. 



□ 



3 The structure function 

Given the random velocity field u{x,-) constructed above, under the assumption of Corol- 
lary|2l a quantity of major interest in the theory of turbulence is the longitudinal structure 
function defined for every integer p and £ > as 

Slie) =me,u{x + ee) -u{x))]P 

where (■, ■) is the Euclidean scalar product in M^, e G is a unit vector and x G 
E, we recall, is the expectation on {Q,A,P). 



(15) 
L^, and 



Remark 5 We warn the reader familiar with the literature on statistical fluid mechanics 
that e here is not the dissipation energy, hut just the spatial scale parameter. In the physical 
literature, it is commonly denoted by i; however, in our mathematical analysis we need 
two parameters: the scale parameter of the statistical observation, which we denote by e, 
and a parameter internal to the model that describes the thickness of the different vortex 
filaments, that we denote by i. 
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The moments Sp{e) depend only on e andp, since u {x, •) is homogeneous and isotropic: 
if e = i? ■ ei where ei is a given unit vector, since the dual of R is R~^, we have 

E [{u{x + ee) - u{x), ef] = E [{u{ee) - m(0), e)^] 

= E[{Ru{eei) -u{Q),eY] 
= E[(M(£ei)-M(0),ein. 

For this reason we do not write explicitly the dependence on x and e. 
Let us also introduce the (non- directional) structure function 

Sp{e) = E[\u{x + ee)-u{x)\''] 

which also depends only on e and p. We obviously have 

Sl{e) < S,{e). 

We shall see that, for even integers p, they have the same scaling properties. At the 
technical level, due to the previous inequality, it will be sufficient to estimate carefully 
sl{e) from below and Sp{e) from above. 



3.1 The main result 

The quantities sl{e) and Sp{e) describe the statistical behavior of the increments of the 
velocity field when e — and have been extensively investigated, see [Hj. They are 
expected to have a characteristic power-like behavior of the form (^, and similarly for 
Sp{e). Our aim is to prove that, for the model described in the previous section with a 
suitable choice of 7, holds true in the sense that the limit 

c^ = ito!2i5M (16) 

^ e^o loge ^ ^ 

exists (similarly for Sp{e)) and is computable. The following theorem gives us the neces- 
sary estimates from above and below, for a rather general measure 7. Then, in the next 
subsection, we make a choice of 7 in order to have the classical multifractal scaling theory. 

Theorem 1 Assume that 

-f{UP£T) < 00 



for every p > 1. Then, for any even integer p > 1 there exist two constants Cp,Cp > 
such that ^ 

cpj [UnrU^e] < she) < Sp{e) < Cp^ ^ 



(17) 



for every e G (0, 1) 
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The proof of this result is long and reported in a separate section. 

We would like to give a very rough heuristic that could explain this result. It must be 
said that we would not believe in this heuristic without the proof, since some steps are 
too vague (we have devised this heuristic only a posteriori). 

What we are going to explain is that 

W(|n,ingie (x + ee) - ^single (x)^ ~ f £T. 



This is the hard part of the estimate. 

Let us discuss separately the case e > t from the opposite one. When e > I the 
vortex structure Msingie is very thin compared to the length e of observation of the dis- 
placement, thus the difference Usingie (x + £e)—Using\e (x) does not really play a role and the 
value of W(|Msingie (x + ee) — Wsingie (x)l^) comes roughly from the separate contributions 
of ^single (x + ee) and Msingie (x), which are similar. Let us compute W(|Msingie (x)|^). 

Consider the expression (jHl) which defines Wsingielx). Strictly speaking, consider the 
short-range case, otherwise there is a correction which makes even more difficult the 
intuition. Very roughly, Ki{x — Xt) behaves like I ■ lxt6B(x/), hence, even more roughly, 
Msingie(x) behaves like 



u r 

'Wsinglc(x) ~ y y '^XteB{x/)dXt. 



When Xt is a smooth curve, say a stright line (at distances compared to £), then ^Xt&B{x/)dXt 
is roughly proportional to i if Xt crosses B {x,i), while it is zero otherwise. We assume 
the same result holds true when Xt is a Brownian motion. In addition, Xt crosses B {x, i) 
with a probability proportional to the volume of the Wiener sausage, which is iT. Sum- 
marizing, we have 

\ i with probability iT 
, ~ I Q otherwise 

Therefore Msingie(x) takes rougly two values, U with probability iT and otherwise. It 

follows that W(|Msingle {x)f) ~ U^T. 

Consider now the case e < i. The difference now is important. Since the gradient of 
Ki is of order one, we have 

U 

■Wsingle (X + ee) - Msinglo (x) ~ J^e J I Xt(iB{x/)dXt. 

As above we conclude that Wgingie (x + ee) — Msingie (x) takes rougly two values, eU / i with 
probability iT and otherwise. It follows that W(|Msingie (x)|^) ~ U'P {jYiT. The intu- 
itive argument is complete. 
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3.2 Example: the multifractal model 

The most elementary idea to introduce a measure 7 on the parameters is to take U and T 
as suitable powers of £, thus prescribing a relation between the thickness i and the length 
and intensity. That is, a relation of the form 

d-f{UJ,T) = 6eH{U)6ia{T)r^di. 

Moreover, we have to prescribe the distribution of i itself, which could again be given by 
a power law i~^di. The K41 scaling described below is such an example. 

Having in mind multi-scale phenomena related to intermittency, we consider a super- 
position of the previous scheme. Take a probability measure 6 on an interval I C M_|_ 
(which measures the relative importance of the scaling exponent h E I). Given two 
functions a,b : I ^ M+ with a{h) < 2 (to ensure i"^ < T) consider the measure 

d-f{U,i,T) = [5^H{U)5^aw{T)r^^^M[\ e{dh). (18) 

Then, according to Theorem [T] , we must evaluate 

7([/Wl£<,) = I £'*P+i+"('^)-*('^)rf£^(d/i) 



while 

7 



[0,£]x/ J[£,+oo)xI 

Cp^h£^P+^^''^''^~''^^^9{dh) 



As e — s> 0, by Laplace method, we get 

lim 1^1^^ = Mp + 3 - D(h)] = C 

with D{h) = h{h)—a{h) + 1. With this choice of 7 we have recovered the scaling properties 
of the multifractal model of [llj. See for a review. 

Consider the specific choice 6{dh) = 6i/3{dh) and a(l/3) = 2, &(l/3) = 4. We have 

Cp=[/ip+3-D(/l)];,=l/3 = |. 

This is the Kolmogorov K41 scaling law for 3d turbulence. The choice a(l/3) = 2, namely 
T = has the following geometrical meaning: the spatial displacement and the thickness 
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of the structure is comparable (remember that the curves are Brownian), hence their shape 
is blob-hke, as in the classical discussions of "eddies" around K41. The choice fe(l/3) = 4, 
namely the measure i~^di for the parameter i, corresponds to the idea that the eddies 
are space-filling: it is easy to see that in a box of unit volume the number of eddies of size 
larger that i is of order Finally, the choice 9(dh) = Si/s{dh), namely U = £^^^, is the 
key point that produces = |; one may attempt to justify it by dimensional analysis or 
other means, but it essentially one of the issues that should require a better foundation. 

Remark 6 K4I can be obtain from this model also for other choices of the functions a{h) 
and h{h). For example, we may take 6{dh) = 6i/3{dh), a(l/3) = and b{l/3) = 2. The 
value choice 6{dh) = 5ii^{dh) is again the essential point. The value a(l/3) = means 
T = 1, hence the model is made of thin filaments of length comparable to the integral 
scale, instead of blob-like objects. The value 6(1/3) = 2, namely the measure i~'^di for 
the parameter I, is again space filling in view of the major length of the single vortex 
structures (recall also that the H aus dor ff dimension of Brownian trajectories is 2). From 
this example we see that in the model described here it is possible to reconcile K4I with a 
geometry made of thin long vortexes. 

Remark 7 Assume inf / > and, for instance, snpjD < 00. Then limp^oo Cp = +00. 
In particular, (p > 3 for some even integer p. Since, by iiy^) and Laplace method, 

for any f3 > 0, the condition of remark^ is satisfied. Therefore the velocity field has a 
continuous modification. 

4 Proof of Theorem U 

Let us introduce some objects related to the structure functions at the level of a single 
vortex filament. Let e be a given unit vector, the first canonical one to fix the ideas. Since 

u{x)= (w^i^gig (x)), then 

{u{ee) - u{0), e) = fi [(^eUsmgie, e)] 

where 

^e'^single •— '^single (^^) ~ ^single (0) 

and similarly 

\u{ee) - U{0) \ = \fi [4Msingle]| • 

Of major technical interest will be the quantities, of structure function type, 

5;(£) = W((4M,i,gie,en 

Sp{e) = W{\6eUsinglef). 

They depend also on i, T, U. 
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4.1 Lower bound 

As a direct consequence of lemma IHl and lemma |2l we have: 
Corollary 3 If k is an odd positive integer, then 



Sl{e)=0. 



Moreover, for any even integer p > 1 there is a constant Cp such that 

S^ie) > CpU [{S ^u, ingle, e)^] 



If we prove that 



Cp7 [S;{e 



for every e G (£, 1) and some constant cj, > 0, then 

7[s;{s)]>c;^[um,^,] 

for every e G (0, 1) and some constant c^' > 0, which implies the lower bound of theorem 
m We have 

S^{e) = j [((^eMsingle, e)^] dXQ. 

Since ^ 
where 

KUy)=Keiy)Ae 
by Burkholder-Davis-Gundy inequality, there is Cp > such that 



where 



dt\Kl{ee - Xt) - K^{0 - Xt 



dxo 



Here K^{y) = {Ke{y),e). Therefore 

TIP r r 

The proof of the theorem is then complete with the following basic estimate. 



^p(') > j 
UP 
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Lemma 8 Given p > 0, there exist Cp > such that for every £^ < T < 1 and 



£ > 



we have 



Proof. Recall that 



Consider the function 



X 



KUx) = -2Q--:r for \x\ >1. 



f {z, a) = Kl{z) - Kl{z + ae) 

= {Kx{z) - Ki{z + ae)) A e 



defined for 2; e and a>l. Let e-*- be any unit vector orthogonal to e. We have 



|/(e^,a)| = 2Q 



= 2g 



^_(e^+ae)Ae 



|e-^ + ae 
1 



le, + ae\ 
>f{e^,l)=CoQ 



for Co = 2 ^1 — 2 2^ . Moreover, Ki is globally Lipschitz, hence there is L > such that 
\f{e^ + w,a)-f{e^,a)\ 

< \Kiie^ + w)- Ki{e^) \ + \Ki{e^ + ae + w)- Ki{e^ + ae) | < L \w\ 

for every w G and a > 1. Therefore, there exists a ball B(e-^,a) C 5(0,2) and a 
constant Ci > such that when z e B{e^, a) we have 

\Kl{z) - Kl{z + ae)\ > Ci 

uniformly in a > 1. Then reintroducing the scaling factor i we obtain that for y e 
B{ee^,ea) 

\Kl{y)-Kl{y + ee)\>C^£ 
uniformly in £ e (0, 1) and e > I. Then we have 
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Hence, if e > £ we have 



I dt\Kl{Q-Xt)-Kl{ee-Xt) 
Jo 



iFTom the lower bound in proved in the next section 

p/2 



The lower bound is proved. 

Remark 8 With a bit more of effort is is also possible to prove the bound 

valid for every e e (0, 1) (not only for e > i). This would be the same as the upper bound, 
but we do not need it to prove that the behaviors as e of the upper and lower bound 
is the same. 



4.2 Upper bound 

Lemma 9 For every even p there exists a constant Cp > such that 

Sp{e) < CpU [\6eUsinglef] 

= Cpi [Sp{e)] . 

Proof. Let [5eUsing\c be the i-th. component of 5eUsing\e (0- We have 



Sp{e)<CpY,E[{li [[5.«single(0]J)1 



2=1 



and thus, by lemma |2l 



i=l 



which implies the claim. 

It is then sufficient to prove the bound 



Sp{e)<CpU^{''-^WT. 



□ 
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Again as above, We have 

•^pi^) = J '^XQ [\5eUsinglef] dXQ 

where, by Burkholder-Davis-Gundy inequahty, there is Cp > such that 

TJP 



where 



Therefore 



We= [ dt\K,i€e-Xt)-K,{0-Xt)\^. 
Jo 



It is then sufficient to prove the bound 



For e > £ it is not necessary to keep into account the closeness of to 0: each term 
in the difference of has aheady the necessary scahng. The hard part of the work has 
been done in lemma El above. 

Lemma 10 Given J9 > 0, there exist Cp > such that for every i"^ < T < 1 and 

e>e 

we have 

Proof. Since 

We < 2W{ee) + 21^(0) 



where W{x) = \K£{x — Xt)f dt, by lemma El we have the result. □ 

. We essentially 

repeat the multi-scale argument in the proof of lemma El with suitable modifications. 



^0 

For e < i need to extract a power of e from the estimate of W 
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Lemma 11 As in the previous lemma, when 

e<i 

we have 

Proof. Since now e is smaller that £ we bound Ki{y) — Kg{z) for |?/ — 2;| < £ as 

\Ke{y) - Ke{z)\ < Ce 

if \y\ < 2i. U \y\ > 2i then > £ and using the explicit form of the kernel we have 
the bound 



\K,{y)-K,{z)\<Cei-] <Ce 



u 



u 



where u is the minimum between \y\ and \z\ and in this case l/u < 1. Then, given a 
partition of [2i, A], say 2i = io < ii < ... < = A, as in the proof of lemmaEl we get 



N 



\Ke{y) - K,{z)f = \Ke{y) - K,{z)f 



<CeH 



i=l 

i=l ^ ' 



l?/l<2^ ^ ^^^-^<\y\<^^ + l|j/|>A 
i=l 

4 / /, \ 4 

2 



£,_i<|y|<£, ^Ce I — ) l|j/|>A 



A 



L|y|>A 



where we have used the fact that {u 

N 



< 2/m for M > 2£. Then 



i=l 



"B(0A)\B(a;,4-i) 



^1 ^B{xAY 



and arguing as in that proof 



A 



2p 



i=l 



p/2 



Then, from lemma El in the form 
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and the obvious bound < T for every Borel set B, we have 



i=l 



and taking the hmit as the partition gets finer: 



W [WP/^] < CeP£P£T + C^M ^ ) ^p^ 



+CeP I (^^) (uA Vt)p{u V Vrfud 



u 



A direct computation of the integral as in the proof of lemma El completes the proof, 
taking the limit as A ^ oo. □ 



5 Auxiliary results on Brownian occupation measure 

In this section we prove the estimates on VV[|L^(.^^.||p/^] which constitute the technical 
core of the previous sections. The literature on Brownian occupation measure is wide, so 
it is possible that results proved here are given somewhere or may be deduced from known 
results. However, we have not found the uniform estimates we need, so we prefer to give 
a full self-contained proofs for completeness. Of course, several ideas we use are inspired 
by the existing literature (in particular, a main source of inspiration has been ^3]). 

First, notice that VV[|Lg(^^ ^-^ 1^/^] does not depend on u. But this is not of great help. 
When p = 2, can be explicitly computed: 



= / d^o I dt dz pt {z - Xq) 

Jr3 Jq J B{u/) 

= dxo dt dzpt{xo) 

JR3 Jo J B{u/) 

= \B{u,£)\ [ dt [ pt{xo)dxo = vji^T 







where is a. geometrical constant and pt {x) is the density of the 3D Brownian motion 
at time t. The estimate of VV[|L^|.^^^|^/^] for general p requires much more work. 
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Let 

TBiu/) = mi{t>0: Xte B{u,i)} 

the entrance time in B{u, £) for tlie canonical process. We continue to denote by Wxq the 
Wiener measure starting at xq and also the mean value with respect to it; similarly for 
W, the (7-finite measure dyVxo{X)dxo. 

Lemma 12 For any p > 0, T > 0, £ > 0, Xo,u e R^, it holds that 



< >V,JrB(,,,) < T]Wo[| 



(20) 



Proof. Let us prove the upper bound. Let us set, for simplicity, r = tb{u,£) A T. 
When T < t < T we have 



Xt e B{u, e)^Xt-Xre B{o, 2e) 



then 



1 



)dt < 



^{{Xt-Xr)eB{0,2e)} 



dt 



^{Xr+t-Xr)eB{0,2e) 



dt 



Then, taking into account that t — T implies = 0, 



^Biu/) < 1{t<T} / l{Xr + t-Xr 

Jo 



)eB{o,2e) 



dt 



which gives us, using the strong Markov property. 



>v.o|i^iM)l"<>v.o 



a 



1t<t / 'i-{x^+t-Xr)eB{o,2e)dt 



Pi 



Wx,[t<T]Wo 



'-XteB{o,2e) 



dt 



-yVxArBiu,e)<T]Wo[{Ll^,,,e)r]- 

The upper bound is proved. 

Let us proceed with the lower bound. Let r' = tb{u,£/2) AT. When t' < t < T we have 



Xt - Xr' e B{0,e/2) ^XtE B{u,e) 
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then 



^ ' '^x^,.^-x^,eB{o,i/2)dt 



^B(u,e) ^ / ^Xt&B{u,e)dt > / lxt-x^,eB{o,e/2)dt 

Jt' Jt' 
i-T-t' 

Jo 

f 

t'<T/2 

Jo 

I ^x^,+t-x^,eB{o,e/2)dt. 
Jo 



rT-T' 

- 

'+* 

fT/2 



,eB{o,e/2)dt 



— lr'<T/2 



Then, using again the strong Markov property with respect to the stopping time r', we 
obtain 

>V.o[l^i(«,^)l1 > mo[rB(u,e/2) < T/2]>Vo[(4g,,/2))1 
and the proof is complete. □ 



Letting p = 1 in the previous lemma and using the scale invariance of BM we obtain 
an upper bounds for yVxo[TB{u,£) ^ T] as 

>V,Jtb(„,,) < T] = >V,„/V2[^B(«/V2,^/V2) < T/2\ < 



and the corresponding lower bound for y^xfi\jB{u,ii2) < T'/S]: 



T 



W.Jr.(„,/, < T/2] = W^.Jr,(^„,/^, < T] > ^-^ot^yV^v/V.)^ 
which leads to the easy corollary: 

Corollary 4 For any p > 0, T > 0, £ > 0, Xq, m G M^, it /io/(i5 that 



^V2xo K(v^«//V2)]"'^o[I^S(0//2) I 



<mo[\Ll^u,e/] 



B{o,V2ey 



Wxo/V2[-^B(n/v/2,v/2£)]^o[|-^i(0,2^) T] 
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Lemma 13 Given a > andp > 0, there exist constants c, C > such that the following 
properties hold true: for every T,i satisfying T/i"^ > a we have 



while ifT/e < a 



Proof. Consider first T/i"^ > a. In distribution 



^Xt&B{0,l)dt 



(22) 



and moreover we have that 



c< Wn 



LXteB(o,i) 



i)dt 



< C 



(23) 



uniformly in T/£^ > a > (the constants depend on a and p). The lower bound is 
obtained by setting T/£^ = a while the upper bound is given by the fact that 



LXte-B(o,i) 



i\dt 



for any p > (see for instance ^3], section 3). From and we get the first claim 
of the lemma. 

Next, if T/(? < a, in distribution: 



^B{0/) — 



^XteB{o,e/VT)d't 



and 



B{o,e/VT)\ 



so the second claim is also proved. 

Lemma 14 Given p > 0, a > 0, there are constants c,C > such that, for T > ai"^ 



while for T < ae^ 



□ 

(24) 
(25) 
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Proof. Let us prove ()24|) . Using lemma IT^ equation ()21|) becomes 



for two constants c, C > (depending on a and p) . Moreover, as we remarked at the 
beginning of the section, 

'B{u/)\ 

Using this identity in (j26|) . we get (f2i|l . 

Now consider eq. ()25p . Assume T < ae^. Using eq. ()2Hl and lemma El we have 

and is a consequence of the identity 

The proof is complete. □ 
The previous lemma solves the main problem posed at the beginning of the section. 
We prove also a related result in finite volume that we need for the complementary results 
of section 12.41 We limit ourselves to the upper bounds, for shortness. 

Lemma 15 Given p > 0, a > 0, n G M? , there is a constant C > and a function 9 (R) 
with lim/j^oo 9 (R) = 0, such that, for R> \u\ + i and T > ai"^ 



L 



\xo\>R \ Vi 



R-{\u\ + e) 



while for T < ae^ 



J\xo\>R \ V 1 / 

Proof. Arguing as in the previous lemma, it is sufficient to prove that 

/ W,, [Li(„,)]rfxo < exp . 
J\xo\>R \ yl J 



^\xo\>R 

We have 



I 'Wxa[L^B(u,t)\dxQ = I dxo dt dzpt{z-xo) 
J\xo\>R ' Jb{o,ry Jo J b{u/) 

3„ dt /" dz 



= zui T J — j I / dxo pt {z - xo) 

< wt^Tg {R- T, u, 



T Jb{u/) \B{u,k,)\ Jb{0,RY 
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with (recall that £ < I) 

g (R; T, u, £) = sup / dxo Pt (z - Xo) 

te{0,T],\z\<\u\+£ J B{0,R)'= 

Now (denoting by (Wt) a 3D Brownian motion) 



g{R;T,u,£)^ sup P{\z-Wt\>R) 
te{o,T],\z\<\u\+e 

< sup P{\Wt\ >R-{\u\+£)) 
te{o,T] 

R-(\u\+£) 

= sup P \Wi\ > ^ 

te{o,T] \ Vt 

<piwi>-^-<i''i^*' 



Vt 

This completes the proof. □ 
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